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Abstract 

A search for pair-produced leptoquarks is performed using e + e~ collision events collected by the 
OPAL detector at LEP at centre-of-mass energies between 189 and 209 GeV. The data sample 
corresponds to a total integrated luminosity of 596 pb _1 . The leptoquarks are assumed to be 
produced via couplings to the photon and the Z°. For a given search channel only leptoquark 
decays involving a single lepton generation are considered. No evidence for leptoquark pair 
production is observed. Lower limits on masses for scalar and vector leptoquarks are calcu- 
lated. The results improve most of the LEP limits derived from previous searches for the pair 
production process by 10-25 GeV, depending on the leptoquark quantum numbers. 
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1 Introduction 



In the Standard Model (SM) quarks and leptons appear as formally independent components. 
However, they show an apparent symmetry in terms of the family and multiplet structure of 
the electroweak interactions. Some theories beyond the SM [1] therefore predict the existence 
of new bosonic fields, called leptoquarks (LQ), mediating interactions between quarks and lep- 
tons. The interactions of leptoquarks with the known particles are usually described by an 
effective Lagrangian that satisfies the requirement of baryon and lepton number conservation 
and respects the SU(3)c <E> SU(2) L <8> U(l)y symmetry of the SM [2,3]. This results in nine 
scalar (S) and nine vector (V) leptoquarks which are colour triplets or antitriplets and are 
grouped into weak isospin triplets (Si and Vi), doublets (5*1/2, S1/2, V1/2 an d V1/2) an d singlets 
(5*0, So, Vq and Vo) 1 . Their properties are shown in Tables 1 and 2. A charge eigenstate within 
a multiplet will be referred to as a "state" and denoted by Si(Q cm ) or Vi(Q em ), where Q em is 
the electric charge in units of e. 

Under these assumptions, only the masses and the couplings to right-handed and left-handed 
leptons, denoted by Ar and Al, remain free parameters, since the couplings to the electroweak 
gauge bosons are completely determined by the electric charge and the third component of the 
weak isospin, while the interactions with gluons are given by the colour charge. Each coupling 
can carry generation indices for the two fermions [3], so that \ tj couples a leptoquark to an i th 
generation lepton and a j th generation quark. In this note only leptoquark decays involving 
a single family of leptons are searched for, while no distinction is made between quarks from 
different generations. This corresponds to the simplifying assumption that A lJ • A mn = if 
i 7^ m. The states with couplings both to right-handed charged leptons and left-handed neutri- 
nos have an unknown branching ratio into a charged lepton and a quark, f3, depending on the 
relative values of the couplings, while for all the other states (5 has a known fixed value. Some 
leptoquarks with couplings to left-handed leptons have the same properties as scalar quarks in 
super symmetric models with R-parity violation [4]. This is the case for S$(— 1/3), Si/2(l/3) 
and S1/2 (—2/3). The results obtained in this analysis can therefore also be interpreted in terms 
of these models. 

Several experimental results constrain theories that predict the existence of leptoquarks. 
Searches for events with leptoquark single production, where a first generation leptoquark 
could be formed as a resonance between an electron 2 and a quark, were performed by the 
ZEUS and HI experiments at the ep collider HERA [5] and by the DELPHI and OPAL exper- 
iments at LEP [6]. Leptoquark masses, M LQ , of 0(100 GeV) are excluded for A values greater 
than O(l0~ 2 ). All LEP and Tevatron experiments have searched for events with leptoquark 
pair production [7-9], setting limits on Mlq as a function of the branching ratio for decay 
into a charged lepton and a quark. The values of these limits range from 99 GeV to 275 GeV 
depending on the decay channel and the spin of the leptoquarks. 



1 In this paper the notation used in [3] is adopted and a scalar multiplet of weak isospin / is denoted Si 
and a vector multiplet Vj. This is slightly different from the notation used in [2] where, on the contrary, the 
multiplcts are denoted by their multiplicity, i.e. S21+1 or V21+1, and different symbols are used for leptoquarks 
with fermion number, F, equal to 2 (S, V) or (R, U). 

2 Charge conjugation is implied throughout this paper for all particles, e.g. positrons are referred to as 
electrons, etc. 
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Table 1: Quantum numbers and couplings for scalar leptoquarks. F = 3B + L is the fermion 
number which is a function of the baryon and lepton numbers B and L, Q ern is the electric 
charge in units of e, J3 is the third component of the weak isospin and (3 is the branching ratio 
of the decay to a charged lepton and a quark of any flavour. Under the assumption of non-zero 
couplings only within a single generation of leptons, u and d denote up- and down-type quarks 
respectively, and the same Table applies to second and third generation leptoquarks with the 
obvious substitutions e — > //, r. In the last column, A 2 is a shorthand for J2j(^ lj ) 2 , where i 
denotes the lepton generation and j the quark flavour. 

In this paper a search is presented for pair-produced scalar and vector leptoquarks of all three 
generations performed with the OPAL detector. Compared to single production by electron- 
quark interactions, pair production has the advantage that all states can be produced, including 
leptoquarks that decay only into a neutrino and a quark. Searches for this channel at LEP are 
able to explore the region of large decay branching ratio into quark-neutrino final states, where 
the Fermilab experiments have reduced sensitivity. The study is based on data recorded during 
the LEP runs from 1998 to 2000 at centre-of-mass energies, y^s, between 189 and 209 GeV. 
The different values of y/s and the corresponding integrated luminosities are listed in Table 3. 

In principle, leptoquarks of all three generations can be pair-produced in e + e~ collisions at 
LEP, by s-channel 7 or Z° exchange and, in the case of first generation leptoquarks, by the 
exchange of a quark in the t- or u-channel [10]. The current upper limits on the couplings A to 
fermions are O(10~ 2 ) in the mass range kinematically accessible; this makes the t- or u-channel 
contribution to the first generation production cross-section (including interference between 
this channel and the s-channel) less than 1% of the pure s-channel contribution. Therefore, in 
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Table 2: Same as Table 1, but for vector leptoquarks. 

the present analysis, only s-channel leptoquark production is considered. Consequently, for a 
given state the cross-section depends on the mass, the electric charge and the third component 
of the weak isospin, but is independent of the A couplings. On the other hand, for couplings 
smaller than (9(10~ 6 ) the lifetime of leptoquarks would be sufficiently long to have interactions 
with the material of detector and to produce a secondary decay vertex, clearly separated from 
the interaction region of the electron beams. This topology is not considered here as the tracks 
of charged particles are required to come from the primary interaction vertex so that, to sum- 
marize, the present analysis covers the region 10~ 6 < A < 1CT 2 . 

The decay of a heavy leptoquark into a charged lepton and a quark leads to final states 
characterized by an isolated energetic charged lepton and a hadronic jet, while for decays into 
a quark and a neutrino, the final state would have large missing energy and a jet. Given the 
assumptions about the A couplings, the following topologies are considered for events resulting 
from the decay of a leptoquark-antileptoquark pair: 

Class A: Two hadronic jets and two neutrinos; it consists of the final states v\ylu^u k and 
v{y[d? d k , where I — e,/i,T and u\ d? are up- and down-type quarks of the j th generation. 
Class B: Two hadronic jets, one neutrino and one charged lepton of the same generation 
(z/ i / ± ^'rf fc ). 

Class C: Two hadronic jets and one pair of oppositely charged leptons of the same generation 
(l + l~uiu k , l+l-d^d k ). 
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YEAR 


(>/*) (GcV) 


fcdt (pb- 1 ) 

J 


1998 


188.6 


169.1 


1999 


191.6 


28.9 




195.5 


72.3 




199.5 


74.7 




201.7 


39.2 


2000 


203.8 


8.5 




205.1 


69.6 




206.3 


63.1 




206.6 


63.8 




208.0 


6.7 


TOTAL 


595.9 



Table 3: Average centre- of- mass energies and corresponding integrated luminosities for the data 
samples used in the analysis. The search for vector leptoquarks includes only the data with 
y/s > 195 GeV. 

2 The OPAL Detector 

The OPAL detector is described in detail in [11]. It was a multi-purpose apparatus having 
nearly complete solid angle coverage 3 . The central detector consisted of a system of tracking 
chambers inside a 0.435 T solenoidal magnetic field as well as of two layers of silicon microstrip 
detectors [12] surrounding the beam-pipe. The tracking chambers included a high-precision 
drift chamber, a large-volume jet chamber and a set of ^-chambers measuring the track co- 
ordinates along the beam direction. The resolution on the transverse momentum of a track 
was given by a pt /p t ~ \J (0.02) 2 + (0.0015 • p t ) 2 / (GeV) 2 and the average angular resolution was 
about 0.3 mrad in <fi and 1 mrad in 9. A lead-glass electromagnetic calorimeter was located 
outside the magnet coil and covered the full azimuthal range for polar angles in the range 
of | cos(#) |< 0.984. It was divided into two regions: the barrel (| cos(#) |< 0.82) and the 
endcaps (| cos(6>)| > 0.81). The energy resolution for high momentum electrons was around 
3%. The magnet return yoke, divided into barrel and endcap sections along with pole tips, 
was instrumented for hadron calorimetry in the region | cos(6>) |< 0.99. Four layers of muon 
chambers covered the outside of the hadron calorimeter. Close to the beam axis the forward 
calorimeter and gamma catcher, together with the silicon-tungsten luminometer [13] and the 
forward scintillating tile counter [14], completed the geometrical acceptance down to 33 mrad 
from the beam direction. 

3 Thc right-handed coordinate system is defined so that the positive direction of the z axis is along the e~ 
beam; r is the coordinate normal to the beam axis, <fi is the azimuthal angle with respect to the positive direction 
of the x-axis (pointing towards the centre of LEP) and 9 is the polar angle with respect to +z. 
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3 Monte Carlo simulations 



At lowest order the s-channel contribution to the differential cross section for the production 
of a pair of scalar leptoquarks of mass M LQ in e + e~ collisions at a centre-of-mass energy yfs is 
given by [10] 

while for vector leptoquarks one has 



a=L,R 



day 
dcos 6 



~8~7 



1 - 4M : 



LQ 



s [ 4 + 



1 _ 3 (l-4M 2 Q /* 



sin 



1 " (l " 4M2 Q / S ) 



a=L,R 



(2) 



where a is the electromagnetic coupling and 

k a (s) = -Q cm + Q z (e) 



s - m 2 + iM z r z 



Q Z (LQ). 



(3) 



Here Q em is the electric charge of the leptoquark, Mz and T% are the mass and the width of 
the Z° boson, and the couplings are given by 



Q Z (LQ) 



Q!(e) = 



Qcm sin 2 6w 



cos 0w sin 9w 
-\ + sin 2 ^ 
cos 9y/ sin 6w ' 



(4) 



Q R (e) = tmidw , 



with I 3 being the third component of the leptoquark weak isospin and 6w the Weinberg angle. 

The Monte Carlo generator LQ2 [15] is used to simulate leptoquark pair events. Initial 
state QED radiation is included. In LQ2 scalar leptoquarks decay isotropically in their rest 
frame, while the angular distribution of decay products of vector leptoquarks depends on the 
helicity state. The hadronization of the final state quark pair is performed by JETSET [16]. 
For scalar leptoquarks, samples of at least 1000 signal events are generated for each value of 
the leptoquark mass from Mlq = 50 GeV to the kinematic limit in steps of 10 GeV or less 
for all the different decay topologies at the centre-of-mass energies with the highest integrated 
luminosities (189, 196, 200 and 206 GeV). The search for vector leptoquarks includes only data 
with y/s > 195 GeV and the signal was simulated for M L q > 70 GeV. Since leptoquarks carry 
colour, they may hadronize before decaying if their couplings to fermions are small. This effect 
is evaluated from Monte Carlo samples of pair-produced scalar quarks decaying via R-parity 
violating couplings. These events have features similar to events of class C and allow the im- 
pact of this effect on the detection efficiencies and on the leptoquark mass reconstruction to be 
estimated and taken into account as a systematic uncertainty. 



All relevant SM background processes are studied using various samples of simulated Monte 
Carlo events for each centre-of-mass energy in the data. Two-fermion events (Z°*/7* — > ff(7), 
with f = q,r and denoted by 2f), are simulated with KK2f [17]. The Monte Carlo programs 
HERWIG [18], PHOJET [19], and BDK [20] are used to generate two-photon (77) events 
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with hadronic and leptonic final states. Other processes with four fermions in the final state, 
4f, including and Z° pair production, are simulated with grc4f [21] and KORALW [22]. 
JETSET [16] is used as the principal model for the hadronization. Besides the main samples 
alternative generators or hadronization models such as KORALZ [23], Vermaseren [24] and 
HERWIG [18], are used to check the expectation from the SM background. Generally, at each 
centre-of-mass energy, the number of simulated events for the background processes corresponds 
to at least fifteen times the integrated luminosity of the data, except for the 77 process where, 
at some centre-of-mass energies, Monte Carlo events corresponding to only about three times 
the data integrated luminosity are available. 

The full response of the OPAL detector [25] is simulated for all the Monte Carlo events. 

4 Analysis 

All the leptoquark event topologies (classes A to C as defined in Section 1) are character- 
ized by large charged track multiplicities and large number of energy deposits (clusters) in the 
calorimeters due to the hadronization of the quark pair. Moreover, in events of classes B and 
C, energetic and well-isolated charged leptons are present. 

The tracks of charged particles reconstructed in the tracking system and the clusters in the 
electromagnetic and hadronic calorimeters are required to satisfy the same quality criteria 
as in [8]. To avoid double counting, calculations of quantities such as visible energy and 
transverse momentum are performed from charged particle tracks and from clusters in the 
electromagnetic and hadron calorimeters following the method explained in [26]. Electron 
and muon identification is performed using standard OPAL algorithms [27]. The electron 
identification is based on the match between the momentum of a track and the energy of a 
cluster in the electromagnetic calorimeter associated to the track; moreover the value of the 
ionization energy loss, AE/dx, measured for the track in the OPAL jet chamber must be in 
agreement with what expected for an electron. The muon identification requires at least two 
hits corresponding to the direction of the track in the muon chambers and minimum energy 
deposition for clusters in the hadron calorimeter associated to the track. The energy of identified 
electrons is given by the energy of the electromagnetic calorimeter cluster, while for muons the 
momentum of the track is used to calculate the energy. Tau lepton identification is performed 
using an artificial neural network algorithm described in detail in [28]. The hadronic jets are 
reconstructed using the Durham algorithm [29]. The resolution on the direction of a jet is about 
25 mrad, while the resolution on the jet energy is 10-20%, depending on the energy itself, the 
jet shape and the detector region. 

4.1 Event preselection 

Several preselection requirements are applied to all classes of events. To reduce the number of 
events due to interactions of the LEP beams with residual gas in the beam-pipe or with its ma- 
terial, at least 20% of the reconstructed tracks are required to satisfy the track quality criteria. 
There must be at least four accepted tracks and at least four accepted electromagnetic clusters 
not associated to any track. Finally, the total visible energy, -E v is, is required to be greater than 
0.25^ and smaller than 1.25-^/s and its fraction deposited in the region | cos(#)| > 0.9 must 
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be less than 50%. 

After the preselection 51218 events are observed in the data and 49690 are expected from 
SM background, mostly from two fermion events. The efficiencies for signal events range from 
86% to 99% for both scalar and vector leptoquarks at all centre-of-mass energies, depending 
on the leptoquark mass and the decay channel. 

4.2 The z/z/qq channel (class A) 

Signal events of class A are characterized by a pair of hadronic jets and large missing energy 
due to the neutrinos escaping detection. The following cuts were applied to the data: 

(Al) The total visible energy, E vis , has to be in the range 0.25 < E vis / \fs < 0.75. 

(A2) Neutrinos or particles escaping along the beam pipe are not detected resulting in a total 
reconstructed momentum vector of the event, p tot , different from the expected value of 0. 
The missing momentum of the event is then defined as p m i SS = — Ptot- The component of 
the missing momentum in the direction transverse to the beam axis, p™ lss , is required to 
be larger than 0.2y/s. 

(A3) Events are required to contain no isolated electron or muon with an energy, E c or E^, 
larger than where the isolation criterion requires that the angle between the 

lepton and the nearest charged track is larger than 10°. The events must also contain no 
tau lepton with an associated output from the neural net used for the identification, O t , 
larger than 0.75. 

(A4) The events are forced into two jets. The angle between the directions of the jets, is 
required to be such that cos(6 l jj) > —0.1. 

(A5) The invariant mass of the two jets, IVLj, has to be smaller than 70 GeV. 

Table 4 shows the numbers of events after each cut, together with the numbers of back- 
ground events predicted from SM Monte Carlo samples, and the efficiencies for signal events 
corresponding to M L q = 90 GeV at = 206 GeV. Cut (A2) greatly reduces the 77 and 2f 
backgrounds. Cuts (A3)— (A5) reject almost completely 77 and 2f events, and are very effi- 
cient against 4f background. In the whole data sample, 28 events survive the selection, while 
22.81^3 (stat.) 4 events are expected from Standard Model processes, with the largest contribu- 
tion, about 40%, due to events with a single W-boson (Wei/). At y/s = 206 GeV the selection 
efficiency for signal events for leptoquarks of mass M L q = 90 GeV is (31.3 ± 0.7(stat.))%. 

Figure 1 shows the distributions of the variables used in the selection for events of class A. 
The discrepancies between the observed data and the expected SM events in the distribution 
of the scaled visible energy E v i s /^/s : Figure 1(a), are related to the bad modelling of Monte 
Carlo 77 events and of the emission of photons in the initial state (initial state radiation). 



4 The statistical error on the expected background is calculated by considering the 68.27% confidence band 
around the number of events surviving the selection, following [30] . 
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Cut 


Data 


Background 


4f 


77 


2f 


e (%) 


(Al) 


28313 


26596.0 


4299.0 


1175.0 


21122.0 


94.1 


(A2) 


1474 


1404.0 


1353.0 


5.5 


45.4 


61.9 


(A3) 


371 


340.1 


313.7 


2.3 


24.1 


52.9 


(A4) 


45 


43.0 


41.3 


0.8 


0.9 


37.3 


(A5) 


28 


22.8 


21.6 


0.8 


0.4 


31.3 



Table 4: The remaining numbers of events after each cut of selection A (Wqq channel) for 
various background processes compared with the observed number of events in the whole data 
sample. The last column contains the signal efficiency for events with Mlq = 90 GeV at 
y/s = 206 GeV. Within the statistical errors the efficiencies for scalar and vector leptoquarks 
are the same, so the mean value is quoted. 



4.3 The Z ± L/qq channel (class B) 

The selection of signal events of class B is different for final states with an electron or muon 
(class Bl) and those with a tau lepton (class B2). 



4.3.1 Electron and muon channels (class Bl) 

(Bl-1) The visible energy must lie in the range 0.5 < E v - ls /^fs < 1.0. 

(Bl-2) The direction of the missing momentum must satisfy | cos(0 miss )| < 0.9. 

(Bl-3) The event is required to contain at least one identified charged lepton (an electron for the 
first generation, a muon for the second). 

(Bl-4) The most energetic charged lepton in the event is considered to be the one produced in 
the decays of the leptoquark pair. The energy and momentum of the escaping neutrino 
are calculated from the missing momentum of the event. The energy of the most energetic 
lepton (the charged lepton or the neutrino) has to be larger than 0.15y/s, while the energy 
of the second one has to be larger than ClOy^. 

(Bl-5) The charged lepton and the neutrino are required to be isolated from other tracks in the 
event by requiring that the angle between each of them and the nearest charged track, 
9 cct or # MjCt for the charged leptons of the first and second generation respectively, #„ )Ct 
for the neutrino, must be at least 10°. 

(Bl-6) The event is forced into two jets after removing the track corresponding to the charged 
lepton. The angle between the jets is required to satisfy cos(#jj) > —0.1. 

(Bl-7) To reject W-pair events, a five constraint kinematic fit is applied, where energy and 
momentum conservation is required and the two-jet system and the two-lepton system 
are constrained to have the same mass, Mjj^t. As the momentum of the neutrino is 
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not measured, the effective number of constraints in the fit is two. Events with a fit 
probability larger than 0.1 and, at the same time, a fitted mass Mjj fit larger than 75 GeV 
are rejected. 

(Bl-8) Finally, to reconstruct the leptoquark mass, a second kinematic fit is applied with the 
same constraints as in cut (Bl-7), but this time pairing the leptons with the jets. Of 
the two possible combinations, the one with the higher fit probability is considered. The 
events are selected if the fitted mass Mlq is larger than 50 GeV and the fit probability 
Pat is larger than 10~ 3 . 

In Table 5 the numbers of events after each cut are shown, together with the numbers of pre- 
dicted background events and the efficiencies for signal events corresponding to Mlq = 90 GeV 
at yfs = 206 GeV. The contribution of 77 events is negligible after cut (Bl-4). Cuts (Bl-4) 
and (Bl-5) are particularly efficient in reducing 2f events. The numbers of events observed in 
the data, 13 for the first generation and 26 for the second, are in agreement with the expecta- 
tion from Standard Model processes, 13.7l^o (stat.) and 24.51^3 (stat.), respectively. About 
80% of the expected background is due to W + W~ events. At y/s = 206 GeV the selection 
efficiencies for signal events for leptoquarks with M L q = 90 GeV are (28.0 ± 1.0(stat.))% and 
(35.8 ± l.l(stat.))% for the first and second generation respectively. 

The distributions of some of the variables used in the selection for class Bl are presented 
in Figure 2 and show a good agreement between the data and the simulated background. 

4.3.2 Tau channel (class B2) 

(B2-1) To account for the additional neutrinos from the tau decay, the total visible energy re- 
quired is smaller than in class Bl: 0.35 < E vis /y/s < 0.85. 

(B2-2) The direction of the missing momentum is required to satisfy | cos(0 miss )| < 0.9. 

(B2-3) The events are required to contain at least one identified tau lepton. 

(B2-4) The tau with the highest value of the output from the neural network algorithm, O r , is 
chosen as the one coming from the decay of the leptoquark pair. Events are accepted if 
O r > 0.75. 

(B2-5) The energy and momentum of the tau candidate are calculated using the tracks associated 
to the tau by the neural network algorithm and all the clusters in the calorimeters within 
a cone of 10° around the track with the largest momentum. Due to the missing energy 
and momentum carried away by the neutrinos produced in the tau decay, the measured 
tau energy cannot be used as an input in kinematic fits. However it is rescaled using 
coefncents obtained by solving the following equation to require energy and momentum 
conservation: 

ClPh + C 2Ph + °3Pr + C4Pmiss = (0, y/s) 
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e ± i/qq 


Cut 


Data 


Background 


4f 


77 


2f 


£ (%) 


(Bl-1) 


37400 


36904.0 


8767.0 


149.0 


27988.0 


94.3 


(Bl-2) 


15413 


15192.0 


6492.0 


11.9 


8688.0 


86.3 


(Bl-3) 


9832 


10152.0 


4270.0 


6.2 


5876.0 


82.9 


(Bl-4) 


2533 


2529.0 


1714.0 


1.9 


813.0 


79.5 


(Bl-5) 


1042 


1121.0 


1109.0 


1.1 


10.4 


70.5 


(Bl-6) 


32 


36.4 


35.9 


0.3 


0.2 


30.0 


(Bl-7) 


17 


18.4 


17.9 


0.3 


0.2 


28.6 


(Bl-8) 


13 


13.7 


13.4 


0.1 


0.2 


28.0 


/x±i/qq 


Cut 


Data 


Background 


4f 


77 


2f 


£(%) 


(Bl-1) 


37400 


36904.0 


8767.0 


149.0 


27988.0 


92.6 


(Bl-2) 


15413 


15192.0 


6492.0 


11.9 


8688.0 


84.3 


(Bl-3) 


5478 


6064.0 


2820.0 


1.0 


3243.0 


80.0 


(Bl-4) 


1311 


1336.0 


1223.0 


< 0.1 


113.1 


76.4 


(Bl-5) 


997 


1044.0 


1040.0 


< 0.1 


3.5 


70.6 


(Bl-6) 


53 


56.2 


55.4 


0.0 


0.8 


37.4 


(Bl-7) 


32 


30.4 


29.6 


0.0 


0.8 


36.7 


(Bl-8) 


26 


24.5 


23.9 


0.0 


0.6 


35.8 



Table 5: Same as Table 4, but for selection Bl (e ± uqq and /i ± z/qq channels). If the number of 
expected events is smaller than 0.1, but still different from 0, the notation "< 0.1" is used. 



where p^ = EjJ, % = 1, 2, are the measured momentum and energy of the jets obtained 
by forcing the event into two jets after having removed all the tracks and clusters belonging 
to the tau, p T = (p T ,E T ) are the same quantities for the tau, and p m i SS = (Pmiss, |Pmiss|) 
is calculated from the missing momentum of the event. The c coefficients are required 
to be positive. The energy of the tau is then taken to be E T ^ t = c 3 E T and events with 
Erfit < ™>t, where m T denotes the nominal mass of the tau lepton, are rejected. The 
momentum of the tau is recalculated using |p Tj fit| = \J {E T ^ t ) 2 — ml and the original 
measured momentum direction. 

The unchanged jet momenta and the rescaled tau momentum are used as inputs to kine- 
matic fits as described in cuts (Bl-7) and (Bl-8) and the events are accepted or rejected 
by the same criteria. Since the energy of the tau is rescaled using energy and momentum 
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conservation, and since the momentum of the neutrino is unmeasured, the fit has only 
one effective constraint. 

After cut (B2-5), cuts similar to (Bl-4)-(Bl-6) are applied by using the energies and momenta 
of the leptons and the jets as obtained from the kinematic fit: 

(B2-6) The energies of the leptons, E T ^ t and for the tau and the neutrino respectively, have 
to satisfy 0.1 < E rfit /^/s < 0.3 and 0.2 < E ufit /^s < 0.4. 

(B2-7) The angle between the tau momentum and the nearest track not belonging to the tau 
candidate is required to be at least 20°. The corresponding angle for the neutrino has to 
be at least 10°. 

(B2-8) The angle between the jets is required to satisfy cos(0jj) > —0.1. 

In Table 6 the numbers of events after each cut are shown, together with the numbers of pre- 
dicted background events and the efficiencies for signal events corresponding to M L q = 90 GeV. 
Cuts (B2-2) and (B2-4) are particularly efficient in rejecting 2f events. 4f events are especially 
reduced by cuts (B2-5) and (B2-8). At the end of the selection 35 events are observed in the 
data, in agreement with the 36.0^^ (stat.) events expected from Standard Model processes, 
about 90% from W^ 1 boson pair production events. At y/s = 206 GeV the selection efficiency 
for signal events for leptoquarks of mass M L q = 90 GeV is (18.0 ± 0.9(stat.))%. 

Figure 3 shows some of the variables used in this selection. The discrepancy between the 
observed data and the simulated SM background in the distribution of the output from the 
neural net, O r , Figure 3(a), is due to an excess of low energy tau candidates in the data. These 
candidates are not selected by cut (B2-6). The excess in the data between 0.3^ and 0.36^ 
for the scaled energy of the tau (Figure 3(c)), predominantly stems from events taken at the 
lowest centre-of-mass energies. 







r± 


uqq 








Cut 


Data 


Background 


4f 


77 


2f 


e(%) 


(B2-1) 


30206 


29212.0 


5884.0 


362.2 


22966.0 


92.5 


(B2-2) 


7049 


6580.0 


4106.0 


57.2 


2417.0 


84.8 


(B2-3) 


6731 


6304.0 


4043.0 


52.5 


2208.0 


83.9 


(B2-4) 


3048 


3038.0 


2803.0 


16.1 


218.5 


64.0 


(B2-5) 


699 


620.4 


561.6 


3.2 


55.6 


40.4 


(B2-6) 


252 


248.0 


236.4 


0.6 


11.0 


26.1 


(B2-7) 


216 


211.1 


206.7 


0.5 


3.9 


24.6 


(B2-8) 


35 


36.0 


34.4 


0.1 


1.5 


18.0 



Table 6: Same as Table 4, but for selection B2 fr ± z/qq channel). 
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4.4 The l + l qq channel (class C) 



Signal events of this type are characterized by the presence of a pair of isolated high energy 
charged leptons of the same generation and opposite charge. The missing energy of the events 
is small for the first and second generation while, in the case of third generation, a significant 
missing energy is expected because of neutrinos produced in the tau decays. Different sets of 
cuts were applied to select events with electrons or muons and to select events with taus. 

4.4.1 Electron and muon channels (class CI) 
(Cl-1) The visible energy is required to satisfy E vis > 0.75-^/s. 

(Cl-2) The presence of at least one pair of identified electrons or muons with opposite charge is 
required. The most energetic leptons of the same generation and of opposite charge are 
called the "pair" in the following. 

(Cl-3) The energy of the most energetic lepton of the pair, E e i or E^i, has to exceed 0.15^/s, 
while an energy E e2 or E^ 2 of at least 0.1^/s is required for the other lepton. 

(Cl-4) An isolation cut is applied by requiring that the angle between each lepton of the pair 
and the nearest charged track is at least 10°. 

(Cl-5) After the exclusion of the tracks corresponding to the lepton pair, the event is forced into 
two jets. Events are rejected if cos(#jj) and cos(#u) are both smaller than —0.8, where 9$ 
and 6*11 are the angles between the two jets and the two leptons respectively. 

(Cl-6) Finally, a kinematic fit with five effective constraints is applied to reconstruct the lepto- 
quark mass by requiring energy and momentum conservation and constraining the two 
lepton-jet pairs to have the same mass. Of the two possible lepton-jet combinations the 
one with the higher fit probability is chosen. Events are accepted if this probability is 
larger than 1CT 6 , while the fitted mass has to be at least 50 GeV. 

The numbers of events after each cut, together with the numbers of expected background events 
and the efficiencies for signal events corresponding to Mlq = 90 GeV, are shown in Table 7. 
Cuts (Cl-2) and (Cl-3) greatly reduce all kinds of background. Cut (Cl-4) totally suppresses 
2f events. The requirements (Cl-5) and (Cl-6) are useful in further reducing four-fermion 
background. In the search for first generation leptoquarks, 20 events are observed in the data 
while 12.81^3 (stat.) are expected from Standard Model background. For the second genera- 
tion, 4 events are observed, the background expectation being 8. 71^7 (stat.). A contribution 
of about 60% to the total background is expected from Z°Z° events. The efficiency of the 
selection for signal events with leptoquarks of mass M L q = 90 GeV is (50.3 ± 0.7(stat.))% and 
(62.8 ± 0.7(stat.))% for the first and second generation, respectively, at y/s = 206 GeV. 

Figure 4 shows some of the variables used to select events belonging to class CI. 

4.4.2 Tau channel (class C2) 

(C2-1) The visible energy must lie in the range 0.45 < E vis /y/s < 0.95. 

(C2-2) The presence of at least one pair of identified taus with opposite charge is required. 
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e+e qq 


Cut 


Data 


Background 


4f 


77 


2f 


£(%) 


(Cl-1) 


22905 


23093.0 


6998.0 


33.0 


16062.0 


95.2 


(Cl-2) 


4687 


5269.0 


1475.0 


7.0 


3787.0 


83.6 


(Cl-3) 


1729 


1969.0 


322.9 


3.1 


1643.0 


79.9 


(Cl-4) 


67 


42.3 


40.4 


0.2 


1.7 


70.4 


(Cl-5) 


50 


32.5 


31.0 


0.1 


1.4 


65.3 


(Cl-6) 


20 


12.8 


12.3 


< 0.1 


0.5 


50.3 


qq 


Cut 


Data 


Background 


4f 




2f 


£(%) 


(Cl-1) 


22905 


23093.0 


6998.0 


33.0 


16062.0 


91.7 


(Cl-2) 


1821 


2047.0 


643.6 


1.0 


1402.0 


80.4 


(Cl-3) 


68 


70.3 


49.3 


0.0 


21.0 


77.4 


(Cl-4) 


29 


28.0 


28.0 


0.0 


0.0 


72.1 


(Cl-5) 


21 


22.6 


22.6 


0.0 


0.0 


67.7 


(Cl-6) 


4 


8.7 


8.7 


0.0 


0.0 


62.8 



Table 7: Same as Table 4, but for selection CI (e + e qq and ji qq channels). 



(C2-3) For each electric charge the tau candidate with the largest output from the neural network 
algorithm is chosen, and the two outputs O ri and O r2 are combined to form the two-tau 
probability : 

try ® r<2 

1 TT - o T1 o T2 +{i-o T1 ){i-o T2 ) 

V TT is required to be at least 0.9. 

(C2-4) As in selection B2, the energy and momentum of each tau of the pair are calculated from 
the tracks associated to the tau by the identification algorithm and all the clusters in 
the calorimeters within a cone of half angle of 10° around the track with the largest mo- 
mentum. The taus are then removed and the event is forced into a two-jet configuration. 
Then an equation similar to the one described in cut (B2-5), but containing the energies 
and momenta of the jets and taus, is solved. 

The unchanged jet momenta and the rescaled tau momenta are used as inputs to the 
kinematic fit described in cut (Cl-6) to reconstruct the leptoquark mass. As the energies 
of the taus are rescaled using energy and momentum conservation, the effective number 
of constraints in the fit is three. The events are selected if the fitted mass is larger than 
50 GeV and the fit probability is larger than 10~ 6 . 
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After cut (C2-4), the following selections similar to (Cl-3)-(Cl-5) are applied using the 
energies and momenta of the taus and the jets obtained from the kinematic fit. 

(C2-5) The energy of the most energetic tau of the pair has to exceed 0.15-y/s, while an energy 
of at least 0.1 y/s is required for the other tau. 

(C2-6) The angle between a tau and the nearest charged track not belonging to the tau itself, 
T i,ct and # T 2,ct, is required to be at least 20° for each candidate. 

(C2-7) Events are rejected if both cos(#jj) and cos(# TT ) are smaller than —0.8. 

The numbers of events after each cut are shown in Table 8, together with the numbers of ex- 
pected background events and the efficiencies for signal events corresponding to M LQ = 90 GeV. 
Cut (C2-3) reduces in particular the background from 2f events. Cut (C2-4) is efficient against 
each kind of background. In the whole data sample 37 events survive the selection, in good 
agreement with the number expected from Standard Model background, of 38.0t|o (stat.), 
mostly due to W ± and Z° pair production processes (about 50% and 20%). At yHs = 206 GeV 
the efficiency for signal events for leptoquarks of mass Mlq = 90 GeV is (33.3 ± 0.7(stat.))%. 

The distributions of some of the variables used in this selection are shown in Fig. 5. 



T+r qq 


Cut 


Data 


Background 


4f 


77 


2f 




(C2-1) 


35243 


34573.0 


7821.0 


187.0 


26565.0 


90.3 


(C2-2) 


20067 


19393.0 


5690.0 


101.3 


13602.0 


81.7 


(C2-3) 


1503 


1506.0 


1290.0 


24.7 


191.4 


61.5 


(C2-4) 


114 


108.1 


94.1 


1.0 


13.0 


39.9 


(C2-5) 


87 


84.7 


75.5 


0.9 


8.3 


39.2 


(C2-6) 


41 


41.1 


37.6 


0.4 


3.1 


34.7 


(C2-7) 


37 


38.0 


34.7 


0.4 


2.9 


33.3 



Table 8: Same as Table 4, but for selection C2 (t + t qq channel). 



5 Results 

No evidence for leptoquark pair production is observed in the data. Table 9 shows the numbers 
of events selected in the data together with the expectations from Monte Carlo simulations of 
the background processes, including the errors, after all cuts for the different signal topologies. 
A separate comparison is made for the data with y/s > 195 GeV because the search for vector 
leptoquarks includes only data collected at these energies. 
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Channel 


189-209 
Data Bkg 


195-209 
Data Bkg 


vvqq 
e ± ^qq 
/ u ± !/qq 

4- 

r^fqq 
e + e~qq 

r+T^qq 


za zz.a_ 13 _ 35 

1 o 1 q 7+2.4+5.8 
IO 10 -' -1.0-5.8 

26 24.5t?; 5 3 tlj 

oi- o ^* r\A-2 7-1-8 1 

35 36.0±f;J±l:i 
20 12.81?;™ 

4 8.7tg;?t£? 
37 38.0^+6.7 


20 15.2l ; 9 l2.3 
10 10.7t ; 7 t 4 4 j 
18 19.6 |;-:; : :; 

ni-r rr-l-1 9-1-6 3 

21 27.7±i;2±6.3 

15 qo + l.7+3.4 
ld y -°-0.9-3.4 

q 7 n +1.6+2.2 
'- u -0.5-2.2 

24 24.6l 2 ; 3 t 4 4 4 4 



Table 9: The numbers of events observed in the data compared to the numbers of expected 
background events from Monte Carlo, in the different search channels considered. The errors 
are statistical and systematic, respectively. 

A total error on the number of expected background events of 19-46% is estimated assuming 
the following sources: 

• The statistical uncertainty due to the limited number of simulated Monte Carlo events 
(8-26%). 

• The uncertainty introduced by the Monte Carlo modelling of the variables used in the 
selections (12-29%). This is evaluated by displacing the cut value on a given variable, x, 
from the original position x$ to a new position x' , to reproduce on the simulated events 
the effect of the cut on the real data. x' is defined by 

X' Q = Oo - (^)data) ^ + (x)bkg 

where (x)d a ta, (^)bkgj ^data and o"bkg are the mean values and the standard deviations 
of the distributions of the variable x for the data and the simulated background. These 
quantities are calculated by the distributions of x given by the events surviving the cuts on 
all the other variables used in the selection. It was checked that using the distributions of x 
at other stages of the selection leads to negligible changes in the values of this uncertainty. 
This procedure is repeated separately for each variable used in the event selections and 
the change in the number of the expected background events due to the displacement 
of the cut is taken as the systematic error from this source. The different contributions 
are added in quadrature. The main contributions are due to the fit probability and the 
reconstructed W boson mass in the selection of events of class B first generation (21%), 
and to the scaled muons' energies in the search for events of class C second generation 
(15%). 

• The error associated with the lepton identification method is evaluated by considering the 
difference between the number of expected events from Monte Carlo background and the 
number of events observed in the data when only the preselection cuts and the request 
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for presence (or absence, for class A) of leptons are made in the different selections. 
Depending on the class of events, this error is found to range from 3% (r ± z/qq channel) 
to 14% (/i ± z/qq channel). 

• Alternative Monte Carlo generators and fragmentation models are used to check the 
number of expected background events. The differences between the numbers obtained 
using these samples and the main Monte Carlo samples are taken as systematic errors 
and are found to contribute a 5-30%, depending on the different selections. 

The error on the integrated luminosity of the data is less than 0.5% at each energy and is 
neglected. 

In Figure 6 the leptoquark mass reconstructed by the kinematic fits is shown for all the 
events surviving the selections for classes B and C, for both the background and a simulated 
signal. For a leptoquark mass Mlq = 90 GeV at the centre-of-mass energy of 206 GeV the 
mass resolution, obtained by a Gaussian fit to the peak region, ranges from 1.3 GeV (/i + /x~qq 
channel) to 5.0 GeV (r ± z/qq channel), while the mean value of the reconstructed mass is be- 
tween 89.8 (r ± i/qq) and 91.6 GeV (e+e'qq). 

The detection efficiencies for the different topologies of signal events, as functions of the 
leptoquark mass Mlq, are listed in Tables 10 and 11 for scalar and vector leptoquarks respec- 
tively, for the centre-of-mass energies where the signal was simulated. 

The systematic uncertainty on the signal efficiency is evaluated to be 8-31% depending on 
the signal topology and the leptoquark mass. This is estimated by taking into account the 
following sources (the quoted errors are relative): 

• The statistical uncertainty due to the limited number of simulated signal events lies in 
the range of 1-28%. 

• In the region between two simulated leptoquark masses, the value of the efficiency is 
calculated by a linear interpolation. The error associated with this procedure is estimated 
to be 2-8%. 

• The uncertainty introduced by the Monte Carlo modelling of the variables used in the 
selections contributes a systematic error between 3 and 29%. The largest relative effects 
are due to the scaled energies of the muons, E^/y/s and E^/x/s, for events of class C, 
second generation, (up to 28%), to the scaled lepton energies, E T ^/\fs and E v fi t /y/s, 
for events of class B, third generation, (up to 25%) and to the fit probability and the 
reconstructed leptoquark mass (up to 19%), for events of the same class. Most of the 
other selection variables contribute uncertainties of less than 10%. 

• The error associated with the lepton identification method (3-14%). 

• The uncertainty due to the flavour of the final state quarks contributes an error of 2- 
8%. This is evaluated by comparing the efficiencies corresponding to all the final states 
with different quark flavours simulated for a given decay channel, characterized by the 
lepton flavour (for example e + e~uu, e + e~dd, e + e~cc, e + e~ss and e + e~bb for class C, first 
generation). The value of the efficiency is taken to be the mean value, and the largest 
difference between the mean and the single contributions is taken as a systematic error. 
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M LQ (GeV) 


50 


60 


70 


75 


80 


85 


90 


95 


99 


102 


Signal topology 


Generation 








s/s = 189 GeV 


Class A 


1,2,3 


11.4 


17.7 


23.9 


26.9 


30.2 


34.4 


37.4 








Class Bl 
Class Bl 
Class B2 


1 

2 
3 


10.0 
13.5 
1.9 


16.1 

22.8 
5.3 


22.7 
29.6 
9.9 


27.3 
33.2 
12.3 


27.9 
35.5 
15.7 


31.6 
39.9 
19.0 


36.0 
41.9 
22.1 








Class CI 
Class CI 
Class C2 


1 

2 
3 


30.2 
35.1 
20.0 


37.4 
42.9 
26.2 


44.5 
52.0 
31.0 


48.2 
57.1 
30.7 


48.3 
61.1 
34.0 


51.6 
65.5 
33.7 


53.7 
66.2 
35.3 












= 196 GeV 


Class A 


1,2,3 


10.3 


16.3 


22.4 


26.4 


28.1 


30.5 


33.8 


37.7 






Class Bl 
Class Bl 
Class B2 


1 

2 
3 


8.9 
11.7 
1.3 


14.7 
21.2 
4.4 


20.9 
28.6 
8.5 


21.9 
31.3 
11.1 


25.6 
34.2 
13.1 


28.6 
36.7 
17.4 


32.3 
41.3 
19.9 


35.3 
43.5 
23.1 






Class CI 
Class CI 
Class C2 


1 

2 
3 


28.9 
33.7 
18.7 


36.1 
41.1 
25.7 


41.8 
49.4 
27.8 


45.1 
52.9 
29.1 


49.7 
61.1 
32.2 


51.2 
62.5 
32.8 


53.3 
65.4 
33.1 


54.3 
67.8 
34.4 










y/s = 200 GeV 


Class A 


1,2,3 


9.7 


15.1 


21.0 


23.7 


26.6 


29.5 


32.1 


34.9 


38.2 




Class Bl 
Class Bl 
Class B2 


1 

2 
3 


8.7 
11.4 
1.2 


14.2 
20.7 
3.9 


18.9 
27.7 
7.1 


23.8 
31.2 
11.3 


25.0 
32.8 
13.0 


29.3 
37.0 
17.7 


30.2 
37.8 
18.3 


34.0 
42.5 
22.9 


35.4 
43.1 
25.4 




Class CI 
Class CI 
Class C2 


1 

2 
3 


27.8 
32.8 
18.0 


35.1 
40.6 
25.0 


39.8 
48.7 
26.9 


44.7 
53.5 
29.4 


46.8 
57.3 
31.7 


48.8 
62.1 
32.0 


51.7 
65.8 
33.3 


54.4 
67.5 
33.6 


55.6 
66.7 
33.6 








y/s = 206 GeV 


Class A 


1,2,3 


9.1 


13.5 


18.9 


21.1 


24.8 


27.4 


30.2 


33.4 


34.2 


38.1 


Class Bl 
Class Bl 
Class B2 


1 

2 
3 


8.6 
10.8 
1.4 


13.8 
20.8 
2.9 


18.2 
25.1 
6.8 


21.2 
28.7 
8.5 


23.7 
32.4 
10.9 


25.8 
34.5 
13.7 


29.1 
36.9 
18.6 


31.9 
38.5 
20.8 


32.9 
39.6 
22.2 


33.1 
41.7 
22.9 


Class CI 
Class CI 
Class C2 


1 

2 
3 


26.0 
31.3 
16.9 


33.8 
40.2 
25.0 


39.5 
45.4 
25.9 


42.6 
51.8 
29.4 


46.8 
57.6 
31.3 


48.0 
59.3 
32.5 


50.0 
63.0 
32.7 


52.9 
66.3 
35.2 


54.1 
65.7 
34.0 


54.4 
66.6 
34.7 



Table 10: The percentage detection efficiencies for scalar leptoquarks for the various selections 
as functions of the leptoquark mass, M L q, and the centre-of-mass energy. 



• In the range of values of the A couplings covered by this analysis the produced lepto- 
quarks may hadronize before decaying. This process is not simulated by the standard 
signal Monte Carlo. The systematic error on the detection efficiencies associated with 
the fragmentation model is estimated to be 2-4%, evaluated by using MC samples with 
pair produced scalar quarks (squarks) with R-parity violating decays. These events have 
features similar to events of class C but in these samples the hadronization step is sim- 
ulated before the squark decay. The efficiencies obtained by applying the selection for 
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M LQ (GeV) 


70 


75 


80 


85 


90 


95 


99 


102 


Signal topology 


Generation 








s/s = 196 GeV 


Class A 


1,2,3 


22.3 


25.5 


28.9 


31.9 


35.9 


40.8 






Class Bl 
Class Bl 
Class B2 


1 

2 
3 


17.9 
25.3 
7.9 


21.3 
28.3 
10.3 


27.0 
32.5 
13.8 


29.5 
35.5 
18.3 


31.2 
38.1 
21.6 


33.4 
41.4 
23.7 






Class CI 
Class CI 
Class C2 


1 

2 
3 


40.5 
49.4 
28.5 


45.5 
55.1 
30.8 


49.2 
59.5 
32.4 


50.1 
62.7 
33.8 


51.4 
64.1 
34.9 


54.5 
66.8 
33.9 










= 200 GeV 


Class A 


1,2,3 


21.4 


24.4 


28.0 


30.8 


34.0 


39.1 


42.7 




Class Bl 
Class Bl 
Class B2 


1 

2 
3 


17.4 
24.5 
7.3 


19.2 
26.9 
9.4 


25.8 
31.4 
12.1 


29.2 
34.8 
17.1 


30.0 
36.7 
20.4 


33.3 
40.6 
23.4 


33.5 
42.3 
23.9 




Class CI 
Class CI 
Class C2 


1 

2 
3 


39.1 
47.7 
27.8 


44.1 
53.8 
30.2 


48.8 
58.1 
32.0 


49.9 
62.4 
33.2 


50.5 
63.3 
35.0 


52.9 
65.6 
34.7 


56.3 
68.0 
33.0 








y/s = 206 GeV 


Class A 


1,2,3 


19.8 


22.0 


25.3 


28.5 


32.4 


34.4 


37.6 


41.1 


Class Bl 
Class Bl 
Class B2 


1 

2 
3 


16.1 
21.5 
6.3 


17.9 
26.5 
8.9 


22.3 
29.6 
11.7 


24.7 
31.3 
14.2 


26.9 
34.7 
17.4 


30.9 
38.9 
21.1 


32.4 
40.4 
22.4 


31.9 
40.7 
24.4 


Class CI 
Class CI 
Class C2 


1 

2 
3 


40.5 
46.2 
27.0 


43.4 
51.0 
28.1 


47.1 
55.6 
31.0 


48.6 
60.0 
34.0 


50.7 
62.6 
33.9 


51.3 
66.6 
34.0 


54.1 
66.4 
35.1 


54.6 
67.2 
34.2 



Table 11: Same as Table 10, but for vector leptoquarks. 



class C to these events are compared to those obtained using the corresponding standard 
leptoquark samples and the differences are taken as the systematic errors. Moreover, for 
the classes of events where the leptoquark mass is reconstructed, the mass distributions 
obtained by using the different samples are also compared and the mean of the absolute 
value of the difference between the contents of corresponding bins in the two distributions 
is taken as a systematic error. This contribution is estimated to be 3-8%. 

• The data sample is divided into 10 energy bins, as shown in Table 3. However the signal is 
not simulated at each energy. At centre-of-mass energies, v^, where no simulation exists, 
the efficiency for a given leptoquark mass is inferred from the sample at the nearest 
simulated energy, \fp . The efficiency is assumed to be the same as the efficiency for the 
mass point at VP with the same Lorentz boost, that is s(y/s, M L q)= s(Vs', \J s' / s M L q). 
The error associated with this assumption is calculated by comparing the efficiencies 
obtained for corresponding masses at the energies at which the signal is simulated. The 
difference is taken as the error and it is found to be 2-7%. 

The polarization of tau leptons from leptoquark decay is not considered in the simulation of 
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tau decay in the signal events belonging to r+r^qq and r z/qq channels. However it has been 
checked that the effect on the detection efficiencies for vector leptoquarks is negligible. All the 
above errors are considered to be independent and added in quadrature. 

For the purpose of setting limits, the events are divided into different search channels by 
considering their centre-of-mass energy, the decay channel and, for events of classes B and 
C, the reconstructed leptoquark mass divided into 1 GeV bins. The confidence level for the 
existence of a signal is calculated following the method described in [31]. A test statistic is 
defined which expresses how signal-like the data are. The confidence levels are computed from 
the value of the test statistic of the observed data and its expected distributions in a large 
number of simulated experiments under two hypotheses: the background-only (b) hypothesis 
and the signal+background (s + b) hypothesis. 

The test statistic chosen is the likelihood ratio, Q, the ratio of the probability of observing 
the data given the s+b hypothesis to the probability of observing the data given the b hypothesis. 
As all the search channels are considered to be statistically independent and to obey Poisson 
statistics, the likelihood ratio can be computed as 

Q=e- stot n(i+s>Ar 

i 

where rii, bi and Sj are the number of observed candidates, the expected background and the 
expected signal in channel % respectively and s tot = J2i s i- 

The confidence level for the b hypothesis is 1 — CL b , representing the fraction of background- 
only experiments which would produce a value of Q more signal-like than the observed data: 

l-CL b = P(Q>Q obs \b). 

If the data agreed perfectly with the expectation from the background-only hypothesis, a value 
of 1 — CL b = 0.5 would be obtained. A lower value indicates an excess of events in the data; a 
higher value indicates a deficit. Similarly, the agreement of the data with the s + b hypothesis 
is tested by the confidence level CL s+b , defined as 

CL s+b = P(Q<Qobs\s + b) 

which can be used to exclude the s + b hypothesis when it has a small value. However, in the 
case of a large downward fluctuation of the observed background, this procedure may exclude 
a signal for which there is no sensitivity. 

To reduce this possibility the ratio 

CL S = CLs+b/CLf, 

is used to set limits instead. A signal is therefore considered excluded at the 95% confidence 
level if CL S < 0.05. 

The expected signal Sj depends on the electroweak quantum numbers of each leptoquark 
and on the unknown leptoquark mass. The assumption is made that for each scenario only one 
state contributes to the cross-section. Therefore, for each state in the model, CL S and 1 — CL b 
must be calculated as a function of Mlq. In the cases of Sq(— 1/3), Si/2(— 2/3), Vq(— 2/3) 
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and Vi/ 2 (— 1/3) the value of the branching ratio into a charged lepton and a quark, (3, is not 
predicted in the model either and exclusion limits are therefore functions of both M L q and (3. 

The statistical and systematic uncertainties on the expected number of events both for sig- 
nal and background are incorporated in the calculations of the confidence levels as suggested 
in [32]. The probability of observing rii events in channel i, and the corresponding value of 
the test-statistic Q, are integrated over possible values of s$ and bi given by their uncertainties, 
assuming Gaussian distributions, with a lower tail cut-off at zero, so that negative Sj or bi are 
not allowed. In this approach the errors on Sj and bi within a channel and between channels 
are considered to be uncorrelated. 

Figures 7-12 and 15-20 show the values of CL S as a function of the leptoquark mass, Mlq, 
for scalar and vector leptoquarks with the branching ratio (3 predicted in the model. The lower 
limit at the 95% CL on Mlq corresponds to the intersection with the line at CL S = 0.05. In 
the same figures the curves representing the values of 1 — CL^ are also shown. In a Gaussian 
approximation a value 1 — CL5 = 4.55 x 10~ 2 would indicate a 2a excess beyond the background 
median expectation and 1 — CL5 = 2.7 x 10~ 3 would indicate a 3<r excess. The vertical scales on 
the right-hand side of Figures 7-12 and 15-20 correspond to this approximation. The regions 
excluded in the (3 — Mlq plane of the states So, Si/2(— 2/3), Vo and V\ji{— 1/3), whose (3 depend 
on the relative weights of the unknown left and right A couplings, are shown in Figures 13, 14, 
21 and 22. The mass limits obtained are summarized in Table 12. Because of the very small 
cross-section and the lower efficiency of the selection for the r ± ^qq channel, this search can 
only improve previous lower limits on the mass of the third generation state So(— 1/3) over a 
part of the (3 range, while for the third generation state Si(— 1/3) with f3 = 0.5 no improvement 
is possible. 

6 Conclusions 

The data collected with the OPAL detector at y/s between 189 and 209 GeV, corresponding to 
a total integrated luminosity of 596 pb -1 , are analysed to search for events with pair produced 
leptoquarks of all three generations. The present analysis covers the region of small values of 
the couplings A to fermions (from 10 -6 to 10~ 2 ). No significant signal-like excess with respect to 
Standard Model predictions is found in the data. Lower mass limits are set for scalar and vector 
leptoquarks under the assumption that, for each scenario, only one leptoquark state contributes 
to the cross-section. The present results improve most of the previous LEP lower limits on 
leptoquark masses derived from searches for events due to the pair production process [7, 8] by 
10-25 GeV, depending on the leptoquark quantum numbers. 



23 



LQ 



So 



So 



Si 



Sl/2 



Sl/2 



V 



V 



V 1/2 



V 1/2 



Qc 



-4/3 



2/3 
-1/3 
-4/3 



-2/3 



-5/3 



1/3 



-2/3 



-5/3 



1/3 
-2/3 
-5/3 



-1/3 



-4/3 



2/3 



-1/3 



-1/3 [0.5,1] 



1 





0.5 
1 



[0,1] 



-2/3 [0.5,1] 



1 





0.5 
1 



[0,1] 



Generation 



69 



79 



450 



99 



100 



98 



69 
100 



97 
79 
101 



45(*) 
99 



94 



100 



94 



100 



93 



98 



89 



97 



99 



96 



99 



99 



97 



102 



102 



101 



99 
102 



101 
99 
102 



97 
101 



99 



102 



99 



102 



98 



101 



101 



99 



101 



99 



Table 12: The 95% CL lower limits on scalar and vector leptoquarks masses, in GeV, as obtained 
from the present analysis. (3 is the branching ratio into a charged lepton and a quark. Limits 
obtained by OPAL using LEP1 data are marked with (*). 
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Figure 1: Class A, the z/z/qq channel: distributions of the selection variables for the data 
(points with bars denoting the statistical error), the estimated Standard Model background 
(filled histogram) and a simulated signal (hatched histogram), with arbitrary normalization, 
corresponding to scalar leptoquarks of mass Mlq = 90 GeV at y/s = 206 GeV. All the distri- 
butions are shown for the events surviving all the cuts applied before the cut on the plotted 
variable, following the description of the selection in the text. The arrows indicate the positions 
of the cuts and the accepted regions. 

(a) The scaled visible energy, (b) The scaled transverse missing momentum, (c) The scaled 
energy of the most energetic lepton (electron or muon), if a lepton is found, (d) The output 
from the neural net, O t , for the tau with the highest value in the event, after cut (A-2). 
(e) The cosine of the angle between the two reconstructed jets, 9^. (f) The invariant mass of 



the two reconstructed jets, M 
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Figure 2: Class Bl, the / ± z/qq channel: same as Figure 1, but for class Bl. All the distributions 
refer to the selection for second generation leptoquarks. 

(a) The absolute value of the cosine of 8 m i SS , the angle between the direction of the missing 
momentum and the z-axis. (b) The scaled energy of the most energetic lepton (fi orv). (c) The 
angle between the direction of the most energetic muon in the event and the nearest charged 
track, (d) The cosine of the angle between the directions of the two reconstructed jets, 9jj. 
(e) The invariant mass, My^t, of the jet-jet system reconstructed by the kinematic fit described 
in cut (Bl-7). The first bin also contains the events failing the fit or with a probability smaller 
than 0.1. (f) The logarithm of the fit probability, Pg t , used to reconstruct the leptoquark mass. 
The first bin also contains the events failing the fit or with a probability smaller than 10~ 15 . 
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Figure 3: Class B2, the r ± z/qq channel: same as Figure 1, but for class B2. 
(a) The output O r from the neural network algorithm, for the tau with the highest output 
in the event, (b) The invariant mass, Mjj^t of the jet-jet system reconstructed by the first 
kinematic fit described in cut (B2-5). The first bin contains also the events failing the fit or 
with a probability smaller than 0.1. (c)— (d) The scaled energies of the tau lepton and the 
neutrino, £ , r ,flt/Vs aD d E u ^ t /y/s, as calculated by the kinematic fit used to reconstruct the 
leptoquark mass, after cut (B2-5). (e)— (f) The angles between the leptons and the nearest 
charged track for the tau, 9 T ,ct, an d the neutrino, 9^ ct , respectively, after cut (B2-6). 
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Figure 4: Class CI , the l + 1 qq channel: same as Figure 1, but for class CI 
refer to the selection for first generation leptoquarks. 

(a) The scaled visible energy, (b) The scaled energy of the second most energetic electron 

in the event, after cut (Cl-2). (c) The angle between the second most energetic electron 

and the nearest charged track, after (Cl-3). (d) The cosine of 9^, the angle between the 

two reconstructed jets, after cut (Cl-4). The distribution does not contain the events with 

cos(6 l ee) > —0.8, which are always selected by cut (Cl-5), independently of the value of cos 

(e) The cosine of 9 ee , the angle between the two most energetic electrons, after cut (Cl-4). 

The distribution does not contain the events with cos(#jj) > —0.8, which are always selected 

by cut (Cl-5), independently of the value of cos(# ee )- (f) The logarithm of the probability of 

the fit used to reconstruct the leptoquark mass. The first bin contains also the events failing 

the fit or with a probability smaller than 10~ 15 . 
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Figure 5: Class C2, the r + r _ qq channel: same as Figure 1, but for class C2. 
(a) The two-tau probability, V TT , defined in cut (C2-3). (b) The logarithm of the probability 
of the kinematic fit used to reconstruct the leptoquark mass. The first bin contains also the 
events failing the fit or with a probability lower than 10~ 15 . (c)-(d) The scaled energies of the 
most energetic and of the second most energetic tau in the event respectively, after cut (C2-4). 
(e)-(f) The angles between the direction of the momenta of the most energetic and the second 
most energetic tau leptons and the nearest charged track, respectively, after cut (C2-5). 
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Figure 6: The leptoquark masses reconstructed by the kinematic fits used in the selections for 
events of classes B and C . The notation is the same as in Figure 1 . 
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CL S Observed CL S Expected for background — — l-CL b Observed 

Figure 7: The value of CL S as a function of the mass, for scalar leptoquarks Si (2/3) and 
Si/2(l/3) with (3 = 0. The observations for the data are shown with solid lines. The shaded 
bands indicate the 68% and 95% probability intervals with respect to the median expectation 
in the absence of a signal (dashed lines). The mass values corresponding to the intersection of 
the observed CL S with the horizontal solid line at CL S = 0.05 represent the exclusion limits at 
95% CL. The dash-dotted line shows the observed values for the Confidence Level 1 — CLf,; its 
median expectation in the background hypothesis (0.5) and the levels for 2a and 3a deviations 
from this value correspond to the horizontal dotted lines. 
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Figure 8: Same as Figure 7, but for Si(— 1/3) with (3 = 0.5. 



33 



OPAL 



— 1 

U 



1 - 



0.05 



10" 2 r 



10" 



t — i — | — i — i — i — i — | — i — i — i — i — | — i — r 

St 



S (-4/3) (3=1, 1" gen. 



i i i i i i i i i ju i i l_ 



t — i — | — i — i — i — i — | — i — i — i — i — | — i — r 
:S (-4/3) (3=1, 2 nd gen. \ 



: i i i i i i i i i i u i i c 



i — i — | — i — i — i — i — | — i — i — i — i — | — i — □ 
:S (-4/3) (3=1, 3 rd gen. 




90 95 100 90 95 100 90 95 100 

M LQ (GeV) M LQ (GeV) M LQ (GeV) 

CL Observed CL Expected for background — — 1-CL. Observed 



Figure 9: Same as Figure 7, but for S (— 4/3) with (3 = 1. 
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Figure 10: Same as Figure 7, but for Si (—4/3) with (3 = 1. 
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Figure 11: Same as Figure 7, but for S%/2(— 5/3) with [3=1. 
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Figure 12: Same as Figure 7, but for Si/ 2 (— 2/3) with (3=1. 
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Figure 13: The 95% CL exclusion curves in the plane (3 vs M LQ for S (— 1/3) with possible 
values of (3 in the range [0.5,1]. The observations for the data are shown with solid lines. 
The shaded bands indicate the 68% and 95% probability intervals with respect to the median 
expectation in the absence of a signal (dashed lines). The excluded region is to the left of the 
solid curve. 
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Figure 14: Same as Figure 13, but for S\/2{— 2/3) with possible values of (3 in the range [0, 1]. 
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Figure 15: Same as Figure 7, but for vector leptoquarks Vi(l/3) and Vi/2(2/3) with [3 = 0. 
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Figure 16: Same as Figure 15, but for V\{— 2/3) with [3 = 0.5. 
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Figure 17: Same as Figure 15, but for V (— 5/3) with {3=1. 
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Figure 18: Same as Figure 15, but for V\{— 5/3) with (3 = 1. 
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Figure 19: Same as Figure 15, but for Vi/2(— 4/3) with (3 = 1. 
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Figure 20: Same as Figure 15, but for Vi/ 2 (— 1/3) with (3 = 1. 
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Figure 21: Same as Figure 13, but for vector leptoquark Vq(— 2/3) with possible values of (5 in 
the range [0.5, 1]. 
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Figure 22: Same as Figure 21, but for V\/2(— 1/3) with possible values of (3 in the range [0, 1]. 
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